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Abstract. The transverse mass spectra of Kg, A, 5 and £1 particles produced in Pb- 
Pb collisions at 40 A GeV/c have been studied for a sample of events corresponding 
to the most central 53% of the inelastic Pb-Pb cross-section. We analyze the 
distributions in the framework of a parameterized model inspired by hydrodynamics. 
The dependence of the freeze-out parameters on particle species and event centrality 
is discussed and comparisons with results at higher energy are shown. 



PACS numbers: 12.38.Mh, 25.75.Nq, 25.75.Ld, 25.75.Dw 
1. Introduction 

Ultra-relativistic collisions between heavy ions are used to study the properties of nuclear 
matter at high energy density. In particular, lattice Quantum Chromo-Dynamics (QCD) 
calculations predict a state of matter of deconfined quarks and gluons (quark-gluon 
plasma, QGP) at an energy density exceeding ~ 0.6 GeV/fm 3 [1]. For recent reviews of 
experimental results and theoretical developments see references [2]. 

Strange particle production has proved to be a powerful tool to study the system 
formed in heavy ion collisions [3]. If a QGP state is formed, an increased production 
of s and s quarks with respect to normal hadronic interactions is expected [4]. The 
abundance of s and s quarks would equilibrate in a few fm/c, comparable with the 
plasma lifetime. The result, after statistical hadronization, would be an enhancement 
of the strange and multi-strange particle abundance in nucleus-nucleus interactions 
with respect to nucleon-nucleon interactions; the enhancement was also predicted to 
increase with the number of valence strange quarks [4]. We have reported recently [5] 
results confirming the WA97 finding [6] that strange particle production in Pb-Pb 
collisions at 158 A GeV/c is enhanced, with the predicted hierarchy as expected in 
a QGP scenario. The effect amounts to about a factor 20 enhancement for the triply- 
strange Q. Preliminary results at 40 A GeV/c [7] indicate a similar pattern of the 
strangeness enhancement to that observed at higher energy. Enhancement of strange 
meson and baryon production in Pb-Pb collisions has also been reported by the NA49 
Collaboration [8, 9, 10]. 

Ratios of particle abundances are rather well described by thermal models [11]; this 
is interpreted as strong evidence of local thermalization of the nuclear matter formed 
in the collision on a time scale shorter than the lifetime of the fireball. Local thermal 
equilibrium is a prerequisite to allow the description of the evolution of the system 
using ideal fluid dynamics, provided that the mean free path of the fireball constituents 
is short enough. Hydro-dynamics or parameterized models inspired by hydro-dynamics 
have shown to give a successful description of a number of observables, i.e. transverse 
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momentum (p T ) and rapidity (y) distributions, direct and elliptic flow, and two-particle 
correlation functions (for recent reviews see, e.g., references [12]). 

Collective dynamics in the transverse direction is of major interest since it can only 
arise by the buildup of a pressure gradient which, in turn, would be strongly suggestive 
of thermal equilibration of the system formed in the collision. 

The shapes of the p T spectra are expected to be determined by an interplay between 
two effects: the thermal motion of the particles in the fireball and a pressure-driven radial 
flow, induced by the fireball expansion. To disentangle the two contributions we rely 
on the blast- wave model [13], which assumes cylindrical symmetry for an expanding 
fireball in local thermal equilibrium, under different hypotheses on the transverse 
flow profile. The analysis of transverse expansion in Pb-Pb at 158 A GeV/c was 
presented in reference [14]. In this paper we discuss the study of the transverse mass 
(m T = \Jp\ + to 2 ) spectra for A, S~, Q~ hyperons, their anti-particles and Kg measured 
in Pb-Pb collisions at 40 A GeV/c. 

2. The NA57 set-up 

The A, S~, Q~ hyperons, their anti-particles and the Kg mesons are identified by 
reconstructing their weak decays into final states containing only charged particles: 

Kg -> 7T+ + 7T- A -> p + n~ 

S" -> A + 7T- n- -> a + k- ^ 

u P + 7T p + 7T 

and the corresponding charge conjugates for anti- hyperons. The charged tracks emerging 
from strange particle decays were reconstructed in a telescope made of an array of 
silicon detector planes of 5x5 cm 2 cross-section placed in an approximately uniform 
magnetic field of 1.4 Tesla perpendicular to the beam line; the bulk of the detectors 
was closely packed in an approximately 30 cm long compact part used for pattern 
recognition. Descriptions of the NA57 apparatus can be found in references [5, 14, 15]. 
There, however, emphasis is placed on the 158 A GeV/c beam momentum set-up. The 
experimental set-up for Pb-Pb collisions at 40 A GeV/c is conceptually similar: the only 
differences are in the telescope inclination and detector positions which have been chosen 
in order to optimize at the lower energy the acceptance and reconstruction efficiency of 
strange particles produced at mid-rapidity. A sketch of the silicon telescope as installed 
on the NA57 optical bench is shown in figure 1. An inclination angle of 72 mrad with 
respect to the beam line and a distance of the first pixel plane from the target equal 
to 40 cm were set in order to accept particles produced in about half a unit of rapidity 
around central rapidity at medium transverse momentum. 

An array of scintillation counters (Petals), placed 7.9 cm downstream of the target, 
provided a fast signal to trigger on the centrality of the collisions. The Petals covered 
the pseudo-rapidity region 0.8 < rj < 1.8 and their thresholds could be set so as to accept 
events with track multiplicities above an adjustable limit. This was tuned so that the 
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Figure 1. Sketch (not to scale) of the NA57 apparatus placed inside the 1.4 T field 
of the GOLIATH magnet in the 40 A GeV/c set-up. All distances are given in mm. 

triggered event sample corresponds to approximately the most central 56% of the Pb-Pb 
inelastic cross-section. The thickness of the Pb target was 1% of an interaction length. 

The centrality of the Pb-Pb collisions is determined (off-line) by analyzing the 
charged particle multiplicity measured by two stations of micro-strip silicon detectors 
(MSD) which sample the pseudo-rapidity intervals 1.9 < rj < 3 and 2.4 < rj < 3.6. 

3. Data sample and analysis 

The results presented in this paper are based on the analysis of the full data sample 
collected in Pb-Pb collisions at 40 A GeV/c, consisting of 240 M events. The selected 
sample used for the analysis corresponds to the most central 53% of the inelastic 
Pb-Pb cross-section. The data sample has been divided into five centrality classes 
(0,1,2,3 and 4, class 4 being the most central) according to the value of the charged 
particle multiplicity measured by the MSD. The procedure for the measurement of the 
multiplicity distribution and the determination of the collision centrality for each class is 
described in reference [16]. The fractions of the inelastic cross-section for the five classes, 
calculated assuming an inelastic Pb-Pb cross-section of 7.26 barn, are the same as those 
defined at higher beam momentum (158 A GeV/c) and they are given in table 1. 

The strange particle signals are extracted with the same procedure used at 158 A 
GeV/c [5, 14] by applying similar geometric and kinematic constraints. The invariant 
mass spectra of 7r~7r + , p7r, Air and AK combinations after all analysis cuts are shown in 
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Table 1. Ccntrality ranges for the five classes. 



Class 


1 2 3 


4 


a/a inel (%) 40 to 53 


23 to 40 11 to 23 4.5 to 11 


to 4.5 



figure 2. The quality of the signals is comparable to that at 158 GeV/c [5, 14], with the 
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Figure 2. Sample invariant mass spectra for tt + tt , p7r, Att and AK. The arrows 
show the PDG mass values [17]. 

hyperon mass peaks at the nominal PDG values [17] and FWHMs of about 15 MeV/c 2 . 
The Kg invariant mass peak has a FWHM of 25 MeV/c 2 and its maximum is shifted 
down by 4 MeV/c 2 with respect to the nominal value. This effect has been accounted 
for in the calculation of the corrections for acceptance and reconstruction inefficiencies. 
The selected particles have been chosen in the invariant mass intervals corresponding to 
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the shaded areas of figure 2. 

The amount of residual combinatorial background has been evaluated using the 
event-mixing technique [14, 18]. As an example, figure 3 shows the it + tt~ invariant mass 
distribution for real and mixed events before (left) and after (right) the application of 
the analysis cuts. For the cascade hyperons, the residual background has been evaluated 




0.4 0.45 0.5 0.55 0.6 0.4 0.45 0.5 0.55 0.6 
Mass( k + 71" ) [GeV] Mass( k + k ) [GeV] 



Figure 3. The 7r + 7r~ invariant mass distribution for a small sample of events before 
analysis cuts (left) and for the total sample after analysis cuts (right). The shaded 
histograms show the combinatorial background evaluated by event mixing. 

to be about 4% for S~ and less than 10% for S + and for the singly-strange 

particles, it has been estimated to be 1%, 0.8% and 2% for Kg, A and A, respectively 

The acceptance regions in the transverse momentum (p T ) versus rapidity (y) plane 
are shown in figure 4. The limits of these windows have been defined in order to exclude 
from the final sample the strange particles whose lines of flight are very close to the 
borders of the telescope, where the systematic errors are more difficult to evaluate. 

Each reconstructed particle is assigned a weight to correct for acceptance and 
reconstruction inefficiencies. The computational algorithm of the event weight is the 
same used at higher energy [5, 14]: a number of Monte Carlo events are generated - 
each event consisting of one simulated particle, with the p T and y of the real particle 
and random azimuthal and internal decay angles, merged with a real event of similar 
telescope hit multiplicity as the original event — and they are reconstructed with the 
same analysis tools as for real events; the final weight is the ratio of the number of 
generated particles to the number of those reconstructed and retained by the selection 
criteria. This individual correction procedure has been applied to all the reconstructed 
Q, H and A particles. For the much more abundant Kg and A candidates, the individual 
weights have been computed for a subsample only. These amount to 1/20 and 1/15 of 
the total sample, respectively, and they have been extracted uniformly over the full data 
taking period. The statistics of particles collected and individually corrected is given in 
table 2. 
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Figure 4. The pi~y acceptance windows. Dashed lines show the position of mid- 
rapidity (y cm = 2.225). 



Table 2. Statistics of the collected and weighted particles. 







A 


A 






total 


97000 


82500 


2100 


439 


39 28 


weighted 
total 


1/20 


1/15 


1 


1 


1 1 



The correction procedure has been extensively checked by comparing real and 
Monte Carlo distributions for several parameters — in particular those used for particle 
selection, such as the distance of closest approach in space between the two decay 
particles, the position of the decay vertices, the impact parameters of the particles^, 
etc. — for different data taking periods and magnetic field orientations (up and down). 

In order to further check the stability of the results (i.e. the — , d2j X distributions) 

J V jtit drm ay > 

\ The impact parameter is approximated as the distance from the primary vertex of the intersection 
of the measured particle trajectory with a plane transverse to the beam line passing through the target 
position. 
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the selection criteria have also been varied, either by changing their limiting values or 
by excluding them one at a time. As a result of these studies we can estimate the 
contribution of the selection and correction procedure to the systematic errors on the 
slope of the distributions to be about 12% for Q and 8% for all other particles. 

The experimental procedure for the determination of the m T distribution is 
described in detail in reference [14], where the results at 158 A GeV/c are discussed. 
The measured distribution of the double differential invariant cross-section — ,^ N , has 

tot arm ay 

been assumed to factorize into a y and an m T dependent part: 

J_J^ toimr) = /to ).J-^ (mr) . (2) 
m T dm T ay m T dm T 

This assumption has been verified by considering the m T (y) distributions for different 
slices in rapidity (transverse mass). 

The shape of the rapidity distribution, i.e. f(y) in equation 2, has been found to 
be well described within our limited acceptance by a constant for all particles except for 
Kg and A, where a Gaussian provides a better description. We have reported a similar 
finding in Pb-Pb collisions at 158 A GeV/c beam momentum [19]. 

The hypotheses on the factorization of the double differential invariant cross-section 
(equation 2) and on the shape of the rapidity distributions (f(y)) can introduce a 
systematic bias on the slope of the — 4^-(m T ) distributions which has been estimated 

J r vm arm y 1 ' 

not to exceed 5%. 

An additional source of systematic errors comes from the residual combinatorial 
background, which can have a different m T distribution than the signal. This error is 
found to be negligible for A and Kg, and about 1.5%, 2%, 4% and 7% for A, S~, S + 
and Q, respectively. 

4. Exponential fits of the transverse mass spectra 

The inverse slope parameter T app ( "apparent temperature" ) has been extracted by means 
of a maximum likelihood fit of the measured double differential invariant cross-section 
— / N A to the formula 

mx ami ay 

1 d*N ( m T \ 

m T dm T dy \ T app J 

The apparent temperature is interpreted as due to the thermal motion coupled with a 
collective transverse flow of the fireball components [13]. The differential invariant cross- 
section distributions are shown in figure 5 as a function of m r for the most central 53% 
of the inelastic Pb-Pb cross-section with the likelihood fit results superimposed. The 
values of the inverse slope parameters T app are given in table 3. The quoted systematic 
errors are obtained by propagating the partial contributions due to i) the selection and 
correction procedure, ii) the hypothesis on the parameterization of the double differential 
cross-section and iii) the residual combinatorial background. 
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Figure 5. Transverse mass spectra of strange particles for the most central 53% of the 
Pb-Pb inelastic cross-section. The superimposed exponential functions have inverse 
slopes equal to the T ap p values obtained from the maximum likelihood fits. 



Table 3. Inverse slope parameter T app (MeV) of the strange particles in the full 
centrality range (0—53%). The first error is statistical, the second one systematic. 







A 


A 






212 ±3 ±21 


261 ±4 ±26 


263 ± 6 ± 26 








n- + n + 


228 ± 12 ± 23 


308 ± 63 ± 31 


368 ± 120 ±40 



The A and A distributions have very similar inverse slopes, agreeing within 2%§. A 
similar baryon-anti-baryon symmetry was reported at top SPS energy for central and 

§ The systematic errors do not play a role in such a comparison since they essentially cancel out in the 
ratio of the spectra. The only different contribution which can be envisaged arises from the different 
residual background. This has been estimated to affect by less than 1.5% the inverse slope of the A, as 
mentioned before. 
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semi-central Pb-Pb collisions for A and S [14]. The similarity of baryon and anti-baryon 
m T slopes is interpreted as suggestive that strange baryons and anti-baryons would be 
produced and evolve in the collision dynamics by similar mechanisms. Based on the A 
and A results presented above, a similar conclusion can be extended to the 40 A GeV/c 
collisions. The large statistical errors associated to the S + m T distribution prevent any 
such conclusion for this multi-strange baryon. 

In the hydro-dynamical view, the apparent temperature depends on m T ||. At a 
given m T0 value, it can be calculated according to the formula [13]: 



^ap P (m T0 ) 



d dN 

lim a — ( ln ;r^) 

TO T ->m T o dm T drrij 



(4) 



With some approximations, this expression simplifies for two asymptotic cases: at 
low p T , i.e. m T — > m , it provides T app = T + |m (/3j_) 2 [13, 20], where T is the 
freeze-out temperature and {/3±) is the average transverse flow velocity; at high p t , i.e. 
mi — > 00, the apparent temperature is simply blue-shifted by the collective dynamics, 
independently of the particle mass [13]: 

Therefore, the apparent temperature actually depends on the m T range where the 
l/m T dN/dm T distribution has been fitted to an exponential function. Nevertheless, 
provided that the spectra are measured at intermediate p T over similar m T — itlq ranges, 
the inverse slopes should follow a hierarchy with the rest mass of the particles. The 
apparent temperature of the E~ hyperon is found to be smaller than those of the A 
and A particles, and of the same order as that of the Kg meson. A similar "violation" 
of the mass hierarchy was reported for the triply-strange f2 particle by the WA97 [21], 
NA57 [14] and NA49 [22] experiments at 158 A GeV/c, and was interpreted [23] as an 
indication of an earlier decoupling of that particle from the expanding fireball. In the 
next section we shall discuss in more details the possibility of a similar effect for the H~ 
at 40 A GeV/c. 

Figure 6 shows a comparison of the inverse slopes measured at the two energies 
plotted as a function of the particle rest mass. The inverse slopes are lower at lower 
energy by about 7% for non-strange and singly-strange particles and about 20% for the 
multi-strange ones. 



4-1. Centrality dependence 

The transverse mass spectra measured in the individual centrality classes of table 1 
are shown in figure 7 for Kg, A, A and S~. For the rarer S + and Q hyperons, the 
collected statistics do not allow to study the centrality dependence. Maximum likelihood 
exponential fits are superimposed to the spectra, as in figure 5. The inverse slope 
parameters T app are given in table 4 as a function of centrality. 

|| In this view, the graphical interpretation of T app would be the inverse of the tangent to the invariant 
1/t7j t d-/V/dra T distribution. See reference [14] for a detailed discussion. 



Study of the m T spectra of strange particles in Pb-Pb at 40 A GeV/c 



11 



1— 1 400 
> 
CD 

^ 375 



CD 350 

Q. 

O 

CO 325 
> 

_E 300 



275 



250 



225 



200 



175 



0.25 0.5 0.75 1 1.25 1.5 1.75 2 2.25 

mass [GeV/c 2 ] 

Figure 6. Apparent temperature as a function of the particle rest mass in Pb-Pb at 40 
A GeV/c (closed symbols, in blue on-line) as compared to 158 A GeV/c (open symbols, 
in red). NA49 results for negative pions [8] and multi-strange hyperons [24, 22, 25] are 
also shown (squares). For display purpose the ordinate scale has been zero-suppressed. 

Table 4. Inverse slopes (MeV) of the to t distributions of Kg, A, A and S~ as a 
function of centrality at 40 A GeV/c. Only statistical errors are shown. Systematic 
errors are estimated to be 10% for all centralitics. 








1 


2 


3 


4 




194 ± 12 


216 ± 7 


209 ±5 


219 ±6 


209 ± 7 


A 


208 ± 12 


255 ±9 


264 ±7 


260 ±7 


279 ±9 


A 


186 ± 17 


248 ± 12 


257 ±11 


303 ± 15 


276 ± 16 




235 ± 67 


230 ± 30 


255 ± 26 


222 ±20 


206 ± 22 



An increase of the apparent temperature when going from class to class 1 is 
observed for the singly-strange particles. Then, from class 1 to 4, the inverse slopes 
of these particles remain approximately constant. The baryon-anti-baryon symmetry 
in the shape of the spectra discussed above for the A hyperon is preserved also as a 
function of centrality. 

No centrality dependence is observed within the errors for the H~ hyperon. 
Remarkably, the apparent temperature of the H~ hyperon is below those of A and 
A and compatible with that of the light Kg meson for all centrality classes apart for 
class where the mass hierarchy appears to be reestablished. 



z 40 158 

r j ^ NA57 

~ ■ § NA49 
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Figure 7. Transverse mass spectra of Kg, A, A and S~ in Pb-Pb collisions at 40 A 
GeV/c for the five centrality classes of table 1. For each species, class 4 is displayed 
uppermost and class lowermost. The spectra of class 2, 3 and 4 have been scaled by 
factor 2, 10 and 100, respectively, for display purposes. 



5. Blast-wave description of the spectra 

In this section we employ the statistical hadronization model of reference [13], which has 
provided a good description of the 158 A GeV/c results [14, 19], to study the strange 
particle m T spectra discussed above. The model assumes that particles decouple from 
a system in local thermal equilibrium with a temperature T, and which expands both 
longitudinally and in the transverse direction; the longitudinal expansion is taken to be 
boost-invariant and the transverse expansion is defined in terms of a transverse velocity 
profile. In this model the differential cross-section for particle j has the form: 

1 ^ - A, j" m lKl ( ^t£) I ( rir (6) 



m T dm T dy J \ T J \ T 
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where p{r) = tanh _1 /^(r) is a transverse boost, K 1 and I are modified Bessel 
functions, Rq is the transverse geometric radius of the source at freeze-out and Aj is a 
normalization constant. The transverse velocity profile f3±(r) has been parameterized 

as 

r < R G (7) 

The numerical value of Rq does not influence the shape of the spectra but just the 
absolute normalization (i.e. the constant Aj). The parameters which can be extracted 
from the fit of equation 6 to the experimental spectra are thus the thermal freeze- 
out temperature T and the surface transverse flow velocity (3s- Assuming a uniform 
particle density, the latter can be related to the average transverse flow velocity using 
the expression (/3±) = 

Equation 6, which is obtained by integrating the Cooper-Frye invariant distribution 
function [26] over the rapidity of fluid elements (77) up to a maximum longitudinal flow 
^max [13], is a good approximation of a full hydro-dynamical calculation for small rapidity 
windows about mid-rapidity. When such a hypothesis is not fulfilled, a numerical 
integration should be performed which, however, requires a priori knowledge of the 
maximum longitudinal flow^f. The smaller the values of the longitudinal flow, the larger 
the deviations from equation 6. In order to estimate the maximum bias introduced by 
this approximation, we have performed the integration by assuming that the longitudinal 
flow has the same strength as the transverse one, whereas the collective expansion is 
expected to be stronger in the longitudinal direction due to the incomplete stopping of 
the incoming nucleons. Both freeze-out parameters T and (/3±) are found to be smaller 
by 2% when the full integration is performed. 

5.1. Global fit with different profiles 

The result of the global fit of equation 6 with a linear profile hypothesis (i.e. with the 
exponent n — I'm equation 7) to the data points of all measured strange particle spectra 
is shown in figure 8 (left panel) for the event sample corresponding to the most central 
53% of the inelastic Pb-Pb cross-section; the fit yields the following values for the two 
parameters T and (fi±): 

T — 118 ± 5(stat)+;j(syst)MeV , {(3 ± ) = 0.392 ± 0.010(stat)+[J;3syst) 

with x 2 /ndf = 77/40. A large contribution to the \ 2 comes from the S spectra: the 
possibility of an early freeze-out of multi-strange particles is discussed below. 

The T and (f3±) parameters are statistically anti-correlated, as can be seen from 
the confidence level contour shown in figure 8 (right panel). The systematic errors on 
T and (fi±) are instead correlated. 

The results of the fits with different profile hypotheses are given in table 5. Contrary 
to the 158 A GeV/c case [14], where the use of the three profiles n=0, 1/2 and 1 results in 

The 7? max can be derived from a fit to the dN/dy distributions as shown, e.g., in reference [19]. 
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Figure 8. Left: blast-wave fits to the transverse mass spectra of strange particles for 
the most central 53% of the Pb-Pb inelastic cross-section at 40 A GeV/c. 
Right: contour plots in the T-(f3±) plane at the ler confidence level as compared to 
the result at 158 A GeV/c [14] in the same centrality range. 

Table 5. Results of the blast- wave model fit using different transverse velocity profiles. 
The quoted errors are statistical. The systematic errors on the temperature and on 
the velocity are estimated to be about 11% and 4%, respectively, for all four profiles. 





n = 


71= 1/2 


n = 1 


n = 2 


T (MeV) 


125 ±4 


121 ±4 


118 ±5 


135 ±10 


<0L> 


0.434 ±0.011 


0.422 ±0.011 


0.392 ±0.010 


0.302 ±0.016 


X 2 /ndf 


41/40 


51/40 


77/40 


167/40 



equally well described spectra" 1 " with similar values of the freeze-out temperature and of 
the average transverse flow velocity, at this energy the quality of the fit is best with n = 
and gradually worsens as the exponent n is increased. Since the case n = is unphysical, 
the choice n = 1/2, which has also been suggested to be a good approximation of the 
full hydro-dynamical calculation [27], appears to be favoured. However, for the sake 
of comparison with the higher energy results, in the following a linear (n — 1) radial 
dependence of the transverse flow velocity is used. 

The la contour plot in the freeze-out parameter space is compared to the 158 
A GeV/c result [14] in the right-hand panel of figure 8. A lower thermal freeze-out 
temperature is measured at the lower beam energy, while the transverse flow velocities 
are found to be compatible within the errors. 

+ The profile n = 2 was instead disfavoured by data. 
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5.2. Particles with/ without quarks in common with the nucleon. 

Since the particles which share valence quarks with the nucleons may have a different 
behaviour from which do not, we have performed separate fits for the particles of the 
two groups. Results of separate blast-wave fits are given in table 6. The freeze-out 



conditions for the two groups of particles are compatible within 2a. Since the interaction 
cross-sections for the particles of the two groups are expected to be very different, this 
finding would suggest limited importance of final state interactions (i.e. a rapid thermal 
freeze-out) and similar production mechanisms for the two groups. A similar conclusion 
has also been drawn at 158 A GeV/c [14], where the freeze-out parameters for the two 
groups agree within la. 

5.3. Earlier freeze-out of multi-strange particles ? 

At higher energies it has been argued that multi-strange baryons may undergo an earlier 
freeze-out than other particles: at top SPS energy this scenario has been suggested 
based on the VL spectra measurements of the WA97 [21], NA57 [14] and NA49 [22] 
Collaborations. At RHIC, results from the STAR Collaboration [28] suggest that the S 
particles freezes out earlier than ir, K, p and A. 

The la contours of the separate blast-wave fits for singly and multiply-strange 
particles are shown in the left-hand panel of figure 9. The plot seems to indicate that 
the S~ hyperon, which statistically dominates the multiply-strange particle fit, may have 
a different thermal freeze-out behaviour from Kg, A and A. However, due to the low 
statistics, it is not possible to extract significant values for both freeze-out parameters 
from the S~ spectrum alone. The possibility of a deviation for the S~ hyperon from 
the freeze-out systematics extracted from the combined fit to the Kg, A and A spectra 
can be better inferred from the integrated information of the S~ spectrum, i.e. from its 
inverse slope. In figure 9 (right panel) we plot a compilation of data* on the m T inverse 
slopes measured in Pb-Pb collisions at 40 A GeV/c, superimposed to the blast-wave 
model results. The full lines represent the inverse slope one would obtain by fitting 
an exponential to a "blast-like" l/mrdN/dm T distribution (i.e. to equation 6) for a 
generic particle of mass m , in the common range 0.05 < m T — m Q < 1.50 GeV/c 2 , for 

* The NA49 results are taken from the following references: ir~ , K + and K~ from [8]; p and deuteron 
from [29]; A and A from [9]; from [30]; S~ from [24] and Q from [22]. The CERES inverse slopes, 
which are taken from reference [31], are shown for negative hadrons, net proton and A. 



Table 6. Thermal freeze-out temperature and average transverse flow velocity in the 
full centrality range. The first error is statistical, the second one systematic. 




K%, A, E~ 120±7± 13 0.39 ±0.02 ±0.02 56.2/21 

A, S + , ir,H + 80 ±22 ±11 0.45 ±0.03 ±0.02 14.8/17 
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two different freeze-out conditions. The first corresponds to the parameters of the best 
fit to the singly-strange particles (T — 115 MeV , (/3±) = 0.40), the second to those 
of the multiply-strange ones (T = 170 MeV, (fi±) = 0.24). Since the inverse slope is a 
function of the m T — m range where the fit is performed, we have also computed the 
blast-wave inverse slopes of Kg, A and S spectra in the m T — mo ranges of NA57 (closed 
circles). The measured values of the inverse slope of the multi-strange baryons appear 
to deviate significantly from the trend of the other strange particles. 

5.4- Centrality dependence 

The hydro-dynamical description of observables related to collective dynamics, e.g. 
the elliptic flow, is strongly influenced by the freeze-out temperature. It is therefore 
important to determine how the thermal freeze-out conditions may change with respect 
to the initial collision geometry. 

The centrality dependence of the freeze-out parameters at 158 A GeV/c beam 
momentum can be summarized as follows: the more central the collisions the larger the 
transverse collective flow and the lower the final thermal freeze-out temperature [14], 
the dependence being more pronounced for the flow. A similar behaviour was reported 
at RHIC by the PHENIX [32], STAR [28] and BRAHMS [33] Collaborations. Higher 
freeze-out temperatures for more peripheral collisions may be interpreted as the result 
of earlier decoupling of the expanding system. 

In figure 10 we show the la confidence level contours for each of the five centrality 
classes defined in table 1. With decreasing collision centrality the transverse flow velocity 




Figure 9. Left: the thermal freeze-out temperature versus the average transverse flow 
velocity for blast- wave fits using a linear [n = 1) velocity profile. The la contours are 
shown, with the markers indicating the optimal fit locations. Right: prediction of the 
blast- wave model for inverse slopes (see text for details). 
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Figure 10. The la confidence level contours from fits in each centrality class. 

decreases steadily, as also observed at higher energies. The freeze-out temperature is 
roughly constant within the errors at the value T « 110 ± 10 MeV for the three most 
central classes (i.e. for the most central 23% of the inleastic Pb-Pb cross-section). 

6. Conclusions 

We have analyzed the transverse mass spectra of A, H~, Q~ hyperons, their anti-particles 
and Kg mesons in Pb-Pb collisions at 40 A GeV/c beam momentum over a centrality 
range corresponding to the most central 53% of the Pb-Pb inelastic cross-section. 

The inverse slopes of these particles are found to be lower than those measured at 
158 A GeV/c beam momentum. As a function of centrality, the inverse slopes of the 
singly-strange particles increase significantly when going from the most peripheral class 
to the next one; for the more central classes, a weak increase with centrality cannot be 
excluded for A and A. The inverse slope of the S~ is compatible with being constant 
over the whole centrality range. S + and Q hyperons have not been studied as a function 
of the centrality due to the limited statistics. 

Particle and anti-particle inverse slopes are compatible within the errors; in 
particular, the inverse slopes of A and A hyperons are the same within 5% over the 
covered centrality range, thus suggesting that strange baryons and anti-baryons may be 
produced and evolve in the collision dynamics by similar mechanisms. 

The inverse slope of the S~ hyperons significantly deviates from the general trend 
of values increasing with the particle rest mass, as observed for non-strange and singly- 
strange particles. 

The analysis of the transverse mass spectra in the framework of the blast-wave 
model suggests that after a central collision the system expands with an average 
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transverse flow velocity of about 40% of the speed of light and then it freezes out when 
the temperature is of the order of 110 MeV. The measured transverse flow velocity is 
compatible with that measured at 158 A GeV/c but the freeze-out temperature is lower 
at low energy. The inverse slope of the S~ hyperon deviates from the value predicted 
by the blast-wave model tuned on singly-strange particles (Kg, A and A). Finally 
the results on the centrality dependence of the expansion dynamics indicate that with 
increasing centrality the transverse flow velocity increases steadily and the freeze-out 
temperature decreases when going from class 1 (23-40%) to the more central classes. 
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